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Abstract 

A search for physics beyond the Standard Model has been performed with high- 
er 2 neutral current deep inelastic scattering events recorded with the ZEUS de- 
tector at HERA. Two data sets, e + p — > e + X and e~p — > e~X, with respective 
integrated luminosities of 112 pb -1 and 16pb _1 , were analyzed. The data reach 
Q 2 values as high as 40000 GeV 2 . No significant deviations from Standard Model 
predictions were observed. Limits were derived on the effective mass scale in 
eeqq contact interactions, the ratio of leptoquark mass to the Yukawa coupling 
for heavy leptoquark models and the mass scale parameter in models with large 
extra dimensions. The limit on the quark charge radius, in the classical form 
factor approximation, is 0.85 ■ 10 -16 cm. 
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1 Introduction 



The HERA ep collider has extended the kinematic range of deep inelastic scattering 
(DIS) measurements by two orders of magnitude in Q 2 , the negative square of the four- 
momentum transfer, compared to fixed-target experiments. At values of Q 2 of about 
4 x 10 4 GeV 2 , the eq interaction, where q is a constituent quark of the proton, is probed 
at distances of ~ 10~ 16 cm. Measurements in this domain allow searches for new physics 
processes with characteristic mass scales in the TeV range. New interactions between e 
and q involving mass scales above the center-of-mass energy can modify the cross sec- 
tion at high Q 2 via virtual effects, resulting in observable deviations from the Standard 
Model (SM) predictions. Many such interactions, such as processes mediated by heavy 
leptoquarks, can be modelled as four-fermion contact interactions. The SM predictions 
for ep scattering in the Q 2 domain of this study result from the evolution of accurate mea- 
surements of the proton structure functions made at lower Q 2 . In this paper, a common 
method is applied to search for four-fermion interactions, for graviton exchange in models 
with large extra dimensions, and for a finite charge radius of the quark. 

In an analysis of 1994-97 e + p data [1], the ZEUS Collaboration set limits on the effective 
mass scale for several parity-conserving compositeness models. Results presented here 
are based on approximately 130 pb _1 of e + p and e~p data collected by ZEUS in the years 
1994-2000. Since this publication also includes the early ZEUS data, the results presented 
here supersede those of the earlier publication [1]. 



2 Standard Model cross section 

The differential SM cross section for neutral current (NC) ep scattering, e^p —>■ e^X, can 
be expressed in terms of the kinematic variables Q 2 , x and y, which are defined by the 
four-momenta of the incoming electron 1 (k), the incoming proton (P), and the scattered 
electron (k') as Q 2 = -q 2 = -{k - k') 2 , x = Q 2 /(2q ■ P), and y = {q • P)/{k ■ P). For 
unpolarized beams, the leading-order electroweak cross sections can be expressed as 



0,Q ' 



dxdQ 2 y ' ^ ' xQ 4 



[l + (l-y) 2 )F^ T (l-(l-y) 2 )xF^] , (1) 



where a is the electromagnetic coupling constant. The contribution of the longitudinal 
structure function, F L (x,Q 2 ), is negligible at high Q 2 and is not taken into account in 
this analysis. At leading order (LO) in QCD, the structure functions F^ 10 and xF^ c are 



1 Unless otherwise specified, 'electron' refers to both positron and electron. 
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given by 



F 2 NC (x,Q 2 ) = Yl [xq(x,Q 2 )+xq(x,Q 2 )] 

q=u,d,s,c,b 

*if c (x,Q 2 ) = Yl B ^ 2 ) [xq(x,Q 2 )-xq(x,Q 2 



q=u,d,s,c 



where q(x, Q ) and q(x, Q ) are the parton densities for quarks and antiquaries. The 



functions A q and B q are defined as 

A(Q 2 ) = \ [(v q L ) 2 + (O 2 + «) 2 + «) 2 ] 

^(Q 2 )= (V q L )(A L q ) - (V q R )(Af) , 
where the coefficient functions V q ,R and A^ ,R are given by: 

V q =Q q - (v e ± a e ) v q xz , 

A* = - (v e ± a e ) a q xz , 
— i 

/ ' 

i q 2 



v f =T^-2sin 2 W Q/ 
«/ = T f 



(2) 



4 sin 2 W cos 2 6 W Q 2 + M 2 Z 



In Eq. (j2J), the superscript z denotes the left (L) or right (R) helicity projection of the 
lepton field; the plus (minus) sign in the definitions of V q and A q is appropriate for 
i = L(R). The coefficients Vf and a/ are the SM vector and axial- vector coupling constants 
of an electron (/ = e) or quark (/ = q); Qf and T'j denote the fermion charge and third 
component of the weak isospin; Mz and 6y/ are the mass of the Z° and the electroweak 
mixing angle, respectively. 



3 Models for new physics 
3.1 General contact interactions 

Four-fermion contact interactions (CI) represent an effective theory, which describes low- 
energy effects due to physics at much higher energy scales. Such models would describe 
the effects of heavy leptoquarks, additional heavy weak bosons, and electron or quark 
compositeness. The CI approach is not renormalizable and is only valid in the low- 
energy limit. As strong limits have already been placed on scalar and tensor contact 
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interactions [2], only vector currents are considered here. They can be represented by 
additional terms in the Standard Model Lagrangian, viz: 

£>ci = Vtjieil^iq^qj) , (3) 

i,j=L,R 
q=u,d,s,c,b 

where the sum runs over electron and quark helicities and quark flavors. The couplings 
rj^j describe the helicity and flavor structure of contact interactions. The CI Lagrangian 
(Eq. (J3J)) results in the following modification of the functions V q l and A l q of Eq. (J2J): 

Q 2 

V q =Qq~ ( V e ± «e) V q Xz + ~Z (ViL + Vm) , 

la 



Q 2 

A \= - ( v e ±a e )a qX z + — {v e d ~ Vm) 



It was assumed that all up-type quarks have the same contact-interaction couplings, and 
a similar assumption was made for down- type quarks 2 : 



Vij Vij Vij j 

n ed = n es = n eb 
'lij 'lij 'lij 7 



leading to eight independent couplings, r/ i , with q — u,d. Due to the impracticality of 



setting limits in an eight-dimensional space, a set of representative scenarios was analyzed. 
Each scenario is defined by a set of eight coefficients, e^J, each of which may take the values 
±1 or zero, and the compositeness scale A. The couplings are then defined by 

Air 

eq _ eq 
lij ij ^2 ' 

Note that models that differ in the overall sign of the coefficients are distinct because 
of the interference with the SM. 

In this paper, different chiral structures of CI are considered, as listed in Tabled Models 
listed in the lower part of the table were previously considered in the published analysis 
of 1994-97 e + p data [1]. They fulfill the relation 

eq . eq eq eq n 

Vll + Vlr - Vrl ~ Vrr = , 

which was imposed to conserve parity, and thereby complement strong limits from atomic 
parity violation (APV) results [3,4]. Since a later APV analysis [5] indicated possible 



2 The results depend very weakly on this assumption since heavy quarks make only a very small con- 
tribution to high-Q 2 cross sections. In most cases, the same mass-scale limits were obtained for CI 
scenarios where only first-generation quarks are considered. The largest difference between the ob- 
tained mass-scale limits is about 2%. 
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deviations from SM predictions, models that violate parity, listed in the upper part of 
Table ^ have also been incorporated in the analysis. The reported 2.3<r deviation [5] 
from the SM was later reduced to around la, after reevaluation of some of the theoretical 
corrections [6, 7]. 

3.2 Leptoquarks 

Leptoquarks (LQ) appear in certain extensions of the SM that connect leptons and quarks; 
they carry both lepton and baryon numbers and have spin or 1. According to the general 
classification proposed by Buchmiiller, Riickl and Wyler [8], there are 14 possible LQ 
states: seven scalar and seven vector 3 . In the limit of heavy LQs [Mlq ^ yfs), the effect 
of s- and t-channel LQ exchange is equivalent to a vector-type eeqq contact interaction 4 . 
The effective contact-interaction couplings, rQ, are proportional to the square of the ratio 
of the leptoquark Yukawa coupling, \lq, to the leptoquark mass, Mlq: 

eq eq ( ^LQ \ ^ 

where the coefficients depend on the LQ species [11] and are twice as large for vector as 
for scalar leptoquarks. Only first-generation leptoquarks are considered in this analysis, 
q = u,d. The coupling structure for different leptoquark species is shown in Table 121 
Leptoquark models Sq and Sy 2 correspond to the squark states dp and ul, in minimal 
super symmetric theories with broken R-parity. 

3.3 Large extra dimensions 

Arkani-Hamed, Dimopoulos and Dvali [12-14] have proposed a model to solve the hierar- 
chy problem, assuming that space-time has 4 + n dimensions. Particles, including strong 
and electroweak bosons, are confined to four dimensions, but gravity can propagate into 
the extra dimensions. The extra n spatial dimensions are compactified with a radius R. 
The Planck scale, Mp ~ 10 19 GeV, in 4 dimensions is an effective scale arising from the 
fundamental Planck scale Mp in D = 4 + n dimensions. The two scales are related by: 

M 2 P ~ R n M 2 D +n . 

For extra dimensions with R ~ 1 mm for n = 2, the scale M& can be of the order of 
TeV. At high energies, the strengths of the gravitational and electroweak interactions can 

3 Leptoquark states are named according to the so-called Aachen notation [9] . 

4 For the invariant mass range accessible at HERA, %fs ~ 300 GeV, heavy LQ approximation is applicable 
for M LQ > 400 GeV. For ZEUS limits covering LQ masses below 400 GeV sec [10]. 
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then become comparable. After summing the effects of graviton excitations in the extra 
dimensions, the graviton-exchange contribution to eq — > eq scattering can be described 
as a contact interaction with an effective coupling strength of [15, 16] 

A 

where Ms is an ultraviolet cutoff scale, expected to be of the order of Md, and the 
coupling A is of order unity. Since the sign of A is not known a priori, both values A = ±1 
are considered in this analysis. However, due to additional energy-scale dependence, 
reflecting the number of accessible graviton excitations, these contact interactions are not 
equivalent to the vector contact interactions of Eq. J3J). To describe the effects of graviton 
exchange, terms arising from pure graviton exchange (G), graviton-photon interference 
(■yG) and graviton- Z [ZG] interference have to be added to the SM eq — > eq scattering 
cross section [17]: 

da{e ± q -> e ± q) da SM | da G | da^ G | da ZG 
dt dt dt dt dt 

= ^^{ 3M4 + 64 ^+ 42 ^ 2+10 ^ 3 + £4 }' 

da^ G ttA aQ q {2u + if 



dt T 2Ml s 2 i 

da ZG ttA a \ (2u + t) 3 i(6u 2 + 6ui + i 2 } 

±v P v„— a P a n 



di 2Mls 2 sm 2 26 w \ e q i — M| eq t - M 2 

where s, t and u, with t = —Q 2 , are the Mandelstam variables, while the other coefficients 
are given in Eq. (J2J. The corresponding cross sections for e ± q scattering are obtained by 
changing the sign of Q q and v q parameters. 

Graviton exchange also contributes to electron-gluon scattering, eg — > eg, which is not 
present at leading order in the SM: 

da(e ± q — > e ± q) ix\ 2 u r A , , 

dt 2M| s 2 1 J 

For a given point in the (x, Q 2 ) plane, the e ± p cross section is then given by 

£a{e 2 d gf X) M 2 ) = ,(,,Q 2 )^ + *,Q 2 )^ + 9(;C ,Q 2 )^, 

where q(x, Q 2 ), q(x, Q 2 ) and g(x, Q 2 ) are the quark, anti-quark and gluon densities in the 
proton, respectively. 
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3.4 Quark form factor 



Quark substructure can be detected by measuring the spatial distribution of the quark 
charge. If Q 2 <C V-^e anc ^ Q 2 ^ V-^V ^ e predictions for the cross sections are 
modified, approximately, to: 

io *»> (. * y / _ 3 < » 



where _R e and are the root-mean-square radii of the electroweak charge of the electron 
and the quark, respectively. 



4 Data samples 

The data used in this analysis were collected with the ZEUS detector at HERA and 
correspond to an integrated luminosity of 48 pb" 1 and 63 pb" 1 for e + p collisions collected 
in 1994-97 and 1999-2000 respectively, and 16 pb _1 for e~p collisions collected in 1998-99. 
The 1994-97 data set was collected at y/a = 300 GeV and the 1998-2000 data sets were 
taken with yfs = 318 GeV. 

The analysis is based upon the final event samples used in previously published cross 
section measurements [18-20]. Only events with Q 2 > 1000 GeV 2 are considered. The 
SM predictions were taken from the simulated event samples used in the cross section 
measurements, where selection cuts and event reconstruction are identical to those ap- 
plied to the data. Neutral current DIS events were simulated using the Heracles [21] 
program with Djangoh [22, 23] for electroweak radiative corrections and higher-order 
matrix elements, and the color-dipole model of Ariadne [24] for the QCD cascade and 
hadronization. The ZEUS detector was simulated using a program based on Geant 
3.13 [25]. The details of the data selection and reconstruction, and the simulation used 
can be found elsewhere [18-20]. 

The distributions of NC DIS events in Q 2 , measured separately for each of the three 
data sets, are in good agreement with SM predictions calculated using the CTEQ5D 
parameterization [26, 27] of the parton distribution functions (PDFs) of the proton. The 
CTEQ5D parameterization is based on a global QCD analysis of the data on high energy 
lepton-hadron and hadron-hadron interactions, including high-Q 2 HI and ZEUS results 
based on the 1994 e + p data. The ZEUS data used in the CTEQ analysis amount to less 
than 3% of the sample considered in this analysis. In general, SM predictions in the Q 2 
range considered here are dominantly determined by fixed-target data at Q 2 < 100 GeV 2 
and x > 0.01 [28]. 
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5 Analysis method 



5.1 Monte Carlo reweighting 

The contact interactions analysis was based on a comparison of the measured Q 2 distri- 
butions with the predictions of the MC simulation. The effects of each CI scenario are 
taken into account by reweighting each MC event of the type ep — > eX with the weight 



d 2 U 

dx dQ' 

W = 



(SM+CI) 



;r%(SM) 

dxdQ 1 \ ) true x . 



(4) 



The weight w was calculated as the ratio of the leading-order 5 cross sections, Eq. ([TJ), 
evaluated at the true values of x and Q 2 as determined from the four-momenta of the 
exchanged boson and the incident particles. In simulated events where a photon with 
energy E 1 is radiated by the incoming electron (initial-state radiation), the electron energy 
is reduced by E^ . This approach guarantees that possible differences between the SM and 
the CI model in event-selection efficiency and migration corrections are properly taken into 
account. Under the assumption that the difference between the SM predictions and those 
of the model including contact interactions is small, higher-order QCD and electroweak 
corrections, including radiative corrections, are also accounted for. 



5.2 Limit-setting procedure 

For each of the models of new physics described above, it is possible to characterize 
the strength of the interaction by a single parameter: 4tt/A 2 for contact interactions; 
{\lq/ M L q) 2 for leptoquarks; A/M| for models with large extra dimensions; and R 2 for 
the quark form factor. In the following, this parameter is denoted by 77. For contact inter- 
actions, models with large extra dimensions and the quark form factor model, scenarios 
with positive and negative 77 values were considered separately. 



For a given model, the likelihood was calculated as 



L( V ) = \e 



where the product runs over all Q 2 bins, ni is the number of events observed in Q 2 bin 
i and ^(77) is the expected number of events in that bin for a coupling strength 77. The 



5 Note that CIs constitute a non-renormalizable effective theory for which higher orders are not well 
defined. 
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likelihood for the complete e^p data set was obtained by multiplying the likelihoods for 
each of the three running periods. 

The value of 77 for which L(rf) is maximized is denoted as rj . First ^ ata ; the value of rj that 
best describes the observed Q 2 spectra was determined. Using ensembles of Monte Carlo 
experiments (MCE), the expected distribution of r] was then determined as a function 
of rjMC, the coupling value used as the input to the simulation. The 95% C.L. limit on rj 
was defined as the value of t]mc for which the probability that \rj \ > |^ ata | was 0.95. 

For each value of t]mci the nominal number of events expected in each Q 2 bin i, denoted 
fiiiVMc) was calculated by reweighting the SM MC prediction according to Eq. (J3J). The- 
oretical and experimental systematic uncertainties were taken into account by treating 
each uncertain quantity as a random variable. For each uncertainty, 100% correlation 
between systematic variations in different bins was assumed. For each individual MCE, 
an independent random variable, 8j, with zero mean, was generated for each systematic 
uncertainty j. The expected number of events in each Q 2 bin i was then given by the 
product of the nominal expectation, /tj, and N sys random factors which account for the 
uncertainties in the estimation of follows: 

N sys 

fJ-i = fi>i{r)Mc) ] I (1 + Cij) Sj . 
3=1 

The coefficent Cjj is the fractional change in the expected number of events in bin i for a 
unit change in 5j . This definition of /i« reduces to a linear dependence of \ii on each 5j when 
5j is small, while avoiding the possibility of \ii becoming negative which would arise if /ij 
was defined as a linear function of the 5j's. For most of the systematic uncertainties, 8j 
follows a Gaussian distribution, except for a few where it follows a uniform distribution, as 
noted in the next section. For a Gaussian Sj distribution, the definition of \ii corresponds 
to a Gaussian distribution in log/ij. About one million MCEs were generated for each 
model, so that the statistical error was negligible. 

5.3 Systematic uncertainties 

Uncertainties in the SM cross sections considered in this study were estimated using the 
Epdflib program [29] based on Qcdnum [30]. Fractional variations estimated from 
Epdflib were used to rescale the nominal SM expectations calculated with CTEQ5D. 
The following uncertainties were included: 

• statistical and systematic uncertainties of the data used as an input to the NLO QCD 
fit. These errors were the largest uncertainty in the SM expectations. At high Q 2 , the 
uncertainty is up to about 4.5% (3%) for e + p (e~p) data; 
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• uncertainty in the value of as(M§) used in the NLO QCD fit. The resulting uncer- 
tainties of NC DIS cross sections at high Q 2 , estimated assuming an error on as(M§) 
of ±0.002 [31], is about 1.6%; 

• uncertainties in the nuclear corrections applied to the deuteron data (Kd) and to 
the data from neutrino scattering on iron (K Fe ) used in Qcdnum. As suggested in 
Epdflib, variations by up to 100% for Kp and 50% for Kp e were applied, treating 
the corrections as uniformly distributed random variables. The corresponding uncer- 
tainties of NC DIS cross sections at high Q 2 , are up to about 1.7% (0.8%) for Kd and 
up to about 3% (0.7%) for K Fe , for e + p (e~p) data. 

The PDF uncertainties calculated using Epdflib are similar to those obtained from a 
ZEUS NLO QCD fit [28], when high-Q 2 HERA data were excluded from the fit. 

In addition to the uncertainty in the SM prediction, the following experimental uncer- 
tainties were taken into account: 

• the scale uncertainty on the energy of the scattered electron of ±(l-3)% depending 
on the topology of the event [32]. The resulting uncertainty of NC DIS cross section 
at high Q 2 is about 0.6% (1.3%), for e + p (e~p) data; 

• the uncertainty in the hadronic energy scale of ±(l-2)% depending on the topology 
of the event [33]. The resulting cross section uncertainty at high Q 2 is about 1%, for 
both e + p and e~p data; 

• uncertainties on the luminosity measurement of 1.6% for the 1994-97 e + p data, 1.8% 
for the 1998-99 e~p data and 2.5% for the 1999-2000 e + p data. Correlations between 
luminosity uncertainties for different data-taking periods are small and were neglected 
in the analysis. 

As the double-angle method used to reconstruct the kinematics of the events [18-20] is 
relatively insensitive to uncertainties in the absolute energy scale of the calorimeter, the 
largest experimental uncertainty in the numbers of NC DIS events expected at high Q 2 
is due to the luminosity measurement. 

6 Results 

No significant deviation of the ZEUS data from the SM prediction using the CTEQ5D 
parameterization of the proton PDF was observed. For all models considered, the best 
description of the data was obtained for very small values of |?7o ata |, i.e. close to the SM. 
The probability of obtaining larger best-fit coupling from the SM, i.e. the probability that 
an experiment would produce a value of \rj \ greater than that obtained from the data, 
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\r] \ > |r/o ata |, calculated with MCEs assuming the SM cross section, was above 25% in all 
cases. Therefore, limits on the strength parameters of the models described in Sec. El are 
presented in this paper. 

The measured Q 2 spectra for e + p and e~p data, normalized to the SM predictions are 
shown in Fig. ^ Also shown are curves, for VV and AA contact-interaction models 
(Section 13. lj) . which correspond to the 95% CL. exclusion limits on A. The 95% C.L. 
limits on the compositeness scale A, for different CI models, are compared in Fig. |2] and 
Tabled Limits range from 1.7 TeV for the LL model to 6.2 TeV for the VV model. Also 
indicated in the figure are the best-fit coupling values, r^ ata = for positive and negative 
couplings. For comparison, the positions of the global likelihood maxima with ±lcr and 
±2cr error 6 bars are included in Fig. El Systematic uncertainties are taken into account 
by averaging the likelihood values over systematic uncertainties. For most models, the 
±2o" error bars are in good agreement with 95% C.L. limits calculated with the MCE 
approach. 

The 95% C.L. lower limits on the compositeness scale A are compared in Table^with lim- 
its from the HI collaboration [34], the Tevatron [35,36] and the LEP [37-40] experiments 
(where only the results from e + e~ — *> qq channel are quoted). In Tabled the relations 
between CI couplings for the compositeness models considered are also included. The 
results on the compositeness scale A presented here are comparable to those obtained by 
other experiments, where they exist. For many models, this analysis sets the only existing 
limits. 

The leptoquark analysis takes into account LQs that couple to the electron and the first- 
generation quarks (u, d) only ( Section l3.2|) . Deviations in the Q 2 distribution of e + p and 
e~p NC DIS events, corresponding to the 95% C.L. exclusion limits for selected scalar and 
vector leptoquark models, are compared with ZEUS data in Fig. El The 95% C.L. limits 
on the ratio of the leptoquark mass to the Yukawa coupling, M L q/\ L q, are summarized 
in Table |2] together with the coefficients a^j describing the CI coupling structure. The 
limits range from 0.27 TeV for S„ mo del to 1.23 TeV for model. Table H also shows 
the LQ limits obtained by the HI collaboration [34] and by the LEP experiments [37,39]. 
In general, comparable limits are obtained. For the S\ , V^j 2 and Vy 2 leptoquarks, the 
ZEUS analysis provides the most stringent limits. 

When only the NC DIS event sample is considered, the leptoquark limits obtained in 
the contact-interaction approximation are similar to, or better than, the high-mass limits 
from the ZEUS resonance-search analysis [10]. However, for Sq , S± and Vq models these 
previously published limits are more stringent, as the possible leptoquark contribution to 

6 Errors are calculated from the likelihood variation: ±ler and ±2cr errors correspond to the decrease of 
the likelihood value to log £(77) = logi(r/ ) — \ and logL(?7) = \o%L{r\a) — 2, respectively. 
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charged current DIS was also taken into account. 

For the model with large extra dimensions ( Section 13. 3|) . 95% CL. lower limits on the 
mass scale in n dimensions of 

M s > 0.78 TeV for A = +1 , 
M s > 0.79 TeV for A = - 1 , 

were obtained. In Fig.HJ effects of graviton exchange on the Q 2 distribution, corresponding 
to these limits, are compared with ZEUS e + p (Fig. 0k) and e~p (Fig. 0b) data. The limits 
on Ms obtained in this analysis are similar to those obtained by the HI collaboration [34] 
and stronger than limits from qq production at LEP [41]. However, if all final states are 
considered, the limits derived from e + e~ collisions exceed 1 TeV [41]. Limits above 1 Te V 
are also obtained in pp from the measurement of e~e + and 77 production [42]. 

Assuming the electron to be point-like (R e = 0), the 95% CL. upper limit on the effective 
quark-charge radius ( Section l3.4|) of 

R g < 0.85 -lO" 16 cm 

was obtained. The present result improves the limits set in ep scattering by the HI 
collaboration [34] (R g < 1.0 ■ 10 -16 cm) and is similar to the limit set by the CDF collabo- 
ration in pp collisions using the Drell-Yan production of e + e~ and /i + yU~ pairs [35] (Rg < 
0.79- 10~ 16 cm). 7 The L3 collaboration has presented a stronger limit (R q < 0.42-10 -16 cm, 
assuming R e = 0), based on quark-pair production measurement at LEP2 [39] and as- 
suming the same effective charge radius for all produced quark flavors. 

If the charge distribution in the quark changes sign as a function of the radius, negative 
values can also be considered for R 2 ,. For such a model, the ZEUS 95% CL. upper limit 
on the effective quark-charge radius squared can be written as: 

-R\ < (1.06- 10~ 16 cm) 2 . 

Cross section deviations corresponding to the 95% CL. exclusion limits for the effective 
radius, R q , of the electroweak charge of the quark are compared with the ZEUS data in 
Fig. 0fc. 

7 Conclusions 

A search for signatures of physics beyond the Standard Model has been performed with 
the e + p and e~p data collected by the ZEUS Collaboration in the years 1994-2000, with 

7 Limits on the effective quark radius published by the CDF collaboration [35] were calculated assuming 
R q = R e . For comparison with limits assuming R e = 0, the limit value was scaled by a factor y/2. 
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integrated luminosities of 112 and 16pb _1 , reaching Q 2 values as high as 4 x 10 4 GeV 2 . 
No significant deviation from Standard Model predictions was observed and 95% C.L. 
limits were obtained for the relevant parameters of the models studied. For the contact- 
interaction models, limits on the effective mass scale, A (i.e. compositeness scale), ranging 
from 1.7 to 6.2 TeV have been obtained. Limits ranging from 0.27 to 1.23 TeV have been 
set for the ratio of the leptoquark mass to the Yukawa coupling, Mlq/Xlq, in the limit of 
large leptoquark masses, Mlq ^> ^/s. Limits were derived on the mass scale parameter in 
models with large extra dimensions: for positive (negative) coupling signs, scales below 
0.78 TeV (0.79 TeV) are excluded. A quark-charge radius larger than 0.85 • 10 _16 cm has 
been excluded, using the classical form-factor approximation. 

The limits derived in this analysis are comparable to the limits obtained by the HI 
collaboration and by the LEP and Tevatron experiments. For many models the analysis 
presented here provides the most stringent limits to date. 
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Table 1: Coupling structure [cll^ £lr: £rl: £rr\ of the compositeness models 
and the 95% C.L. limits on the compositeness scale, A, resulting from the ZEUS 
analysis of 1994-2000 e ± p data. Each row of the table represents two scenarios 
corresponding to t] > (A + ) and rj < (A~ ). The same coupling structure applies 
to d and u quarks, except for the models Ul to U6, for which the couplings for 
the d quarks are zero. Also shown are results obtained by the HI collaboration, the 
pp collider experiments D0 and CDF, and the LEP experiments ALEPH, L3 and 
OPAL. For the LEP experiments, limits derived from the channel e + e~ —> qq are 
quoted. 
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Table 2: Coefficients a?* defining the effective leptoquark couplings in the contact- 
interaction limit Mlq ^> -y/i and the 95% C.L. lower limits on the leptoquark mass 
to the Yukawa coupling ratio M^q/X^q resulting from the CI analysis of the ZEUS 
1994-2000 e ± p data, for different models of scalar (upper part of the table) and 
vector (lower part) leptoquarks. Also shown are results obtained by the HI collabo- 
ration and corresponding contact-interaction limits from the LEP experiments L3 
and OPAL. The limits from LEP on the compositeness scale A, for models with 
coupling structure corresponding to those of scalar (vector) leptoquarks, were scaled 
by factor 1 / \/8tt (1/VAn). 
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Figure 1: ZEUS data compared with 95% C.L. exclusion limits for the effective 
mass scale in the VV and A A contact-interaction models, for positive (A. + ) and 
negative (A~ ) couplings. Results are normalized to the Standard Model expectations 
calculated using the CTEQ5D parton distributions. The insets show the comparison 
in the Q 2 < 10 4 GeV 2 region, with a linear ordinate scale. 
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Figure 2: Confidence intervals o/±l/A 2 at 95% C.L. /or general CI scenarios 
studied in this paper (dark horizontal bars). The numbers at the right (left) margin 
are the corresponding lower limits on the mass scale A + (A~). The dark filled (open) 
circles indicate the positions corresponding to the best-fit coupling values, ^ ata ; for 
positive (negative) couplings. The light filled circles with error bars indicate the 
position of the global likelihood maximum. For calculation of ±ler and ±2<r errors 
on the global maximum position, likelihood values are averaged over systematic 
uncertainties. 
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Figure 3: ZEUS data compared with 95% C.L. exclusion limits for the ratio of 
the leptoquark mass to the Yukawa coupling, M/X, for the Sy 2 , Sf, Vf and Vjy 2 
leptoquarks. Results are normalized to the Standard Model expectations calculated 
using the CTEQ5D parton distributions. The insets show the comparison in the 
Q 2 < 10 4 GeV 2 region, with a linear ordinate scale. 
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Figure 4: ZEUS e + p data (a) and e~p data (b) compared with 95% C.L. exclusion 
limits for the effective Planck mass scale in models with large extra dimensions, 
for positive (M£) and negative (M$) couplings, (c) Combined 1994-2000 data 
compared with 95% C.L. exclusion limits for the effective mean-square radius of 
the electroweak charge of the quark. Results are normalized to the Standard Model 
expectations calculated using the CTEQ5D parton distributions. The insets show 
the comparison in the Q 2 < 10 4 GeV 2 region, with a linear ordinate scale. 
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